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Abstract 

We explicitly show that, in a system with T-duality symmetry, the conhg- 
uration space volume degrees of freedom may hide on the surface boundary 
of the region of accessible states with energy lower than a hxed value. This 
means that, when taking the decompactihcation limit (big volume limit), a 
number of accessible states proportional to the volume is recovered even if no 
volume dependence appears when energy is high enough. All this behavior 
is contained in the exact way of computing sums by making integrals. We 
will also show how the decompactihcation limit for the gas of strings can be 
dehned from a microcanonical description at hnite volume. 
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1 Introduction 


The density of states or its integrated version measuring the number of ac¬ 
cessible sates for a single object (we will denote it as T and reserve hi for the 
density of states) can be considered as the building block to get a description 
of a system with a big number N (of the order of Avogadro’s number) of such 
objects. It is precisely this quantity which provides us with the criterion to 
know whether classical Maxwell-Boltzmann (MB) counting can be applied 
to a system by inspecting whether the number of accessible states for one 
object is much bigger than the number of objects in the system. 

With the density of states for a single object, it is immediate to compute 
q = exp {—Pei) by converting the sum into one over energy levels. 

With q, the partition function for the system with N objects can be calculated 
by using quantum or classical (MB) statistics. 

To compute T{E, V, 1), where we explicitly note that the number of ob¬ 
jects is one, it is necessary to previously look microscopically at the problem 
of the single object at hand in some conhguration in which a hnite volume can 
be dehned. In general, this volume will be big as compared with some other 
length or with energy if we have a fundamental constant with dimensions of 
length at hand. Traditionally, dealing with objects that can be abstracted 
as point particles, this is achieved by studying the quantum mechanical dy¬ 
namics of the object in a box of volume V = where the walls of the box 
are represented by inhnite barriers of potential. However our problem is one 
in which the generic single object is a string. We take the string as funda¬ 
mental in various senses, but particularly in a sense in which everything is 
made of or comes from strings. There is an immediate definition of volume if 
we take our string and put it by hand in a space which is compact ab initio, 
for example, if the space is a hypertorus of volume V = {2tiR)^ where R 
is the value of the radius for every spatial dimension as a simple election. 
Volume is then dehned as a property of the target space in which the string 
propagates. This is the kind of box we are using here. 

Strings are extended objects. As such, they have some kind of non local 
features. What happens to them in compact spaces when holes are present 
is that winding is possible. Winding allows strings to know dimensions as 
a whole. But, quantum mechanically, the wave function has also a kind of 
non-local character that allows even a free particle to know that there is a 
circular dimension and there is also a winding here measuring the number 
of times the center of mass (one point) winds around the circle. The string 
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can simultaneously wind in the two fashions so that we cannot distinguish 
one kind of winding from the other and finally we cannot discern a big circle 
from a small one because each type of winding respectively contributes to the 
energy in a way which is directly and inversely proportional to the radius. 

Historically, the computations of V, 1) have been indirect by using 
canonical (really macrocanonical with chemical potential p = 0) calculations 
to get this quantity from the partition function for a single string by inverse 
Laplace transformation on the inverse temperature variable. Final results 
have been gotten as an asymptotic approximation valid for high energy (cf., 
for instance, P] and [2]). To compute Q{E, V, 1) one has to sum over winding 
and momentum number, when the approximation of substituting the sum 
over them by a multiple integral holds, the output is a series over the energy 
independent of the volume. From the sum itself it is also easy to see that 
there is a regime which is proportional to the volume. Then, two important 
questions emerge: for the single density of states, can we connect the high 
energy regime independent of the volume and the one which is proportional 
to it? For the gas of many strings, is there a well dehned relation among 
energy and volume for which one can say that the system behaves as the gas 
of strings in uncompactified space? We will see that these two questions are 
actually entangled and closely related to the way in which the system passes 
from the low energy regime to the high energy situation as explained in j2]. 

The starting point of the analysis will be the simple form of F as a count¬ 
ing of states with energy lower than a prescribed generic value. r{E,R, 1) 
will be computed in section 2. There, we will explain the easy way in which 
the decompactihcation limit can be dehned and obtained from a simple re¬ 
lation between R and E and how it can be connected to the description of 
the well known particle gas. We will also explain the difficulties of applying 
an analogous reasoning to try to get the dependence on R (or the volume) 
from the sum over windings and momenta when there is no such simple rela¬ 
tion between energy and volume. Furthermore, we will also justify the need 
to quantitatively know more about the approximation of sums by integrals 
thorough the use of the Euler-Maclaurin (EM) formula. So, section 3 we will 
be devoted to explain what is relevant about that celebrated formula. In 
section 4, there appear applications with special emphasis in showing up the 
way the volume independent situation can be corrected thanks to the EM 
formula in order to hnally get the regime in which r{E,R,l) is proportional 
to the volume. In section 5, we will present and explain the dehnition of the 
decompactihcation regime for the gas with an Avogadro’s number of strings 
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clearly showing that a decompactihcation limit can be dehned upon the mi- 
crocanonical description of the system. Finally, section 6 will present some 
comments and possible future prospects for the subject. 

2 Exact forms of P 

Our departure point will be the formulae giving the characteristic perturba¬ 
tive spectrum of a closed string in a toroidal target space, namely 

a'+ 2N + 2N (2.1) 

a' 

N — NEfh-n = ^ (2.2) 

As usual, N and N stand for the right and left moving oscillator numbers and 
start from —1 if we deal with the closed bosonic string and from zero for the 
closed superstring. So m and n are twenty hve dimensional integer vectors 
in the bosonic case and nine dimensional integer ones when supersymmetry 
is present. The second equation is the left-right level matching condition, 
which now, in the compact case, whenever hnite winding and momentum are 
present, allows different right and left moving oscillator numbers. 

With this, it is easy to write T{E,R, 1) which represents the number of 
states accessible to one string with energy lower than E. In a hypertorus 
with d (the number of spatial dimensions, i.e. d = 9 or 25) circles with a 
common radius R related to the volume hj R = I/^/'^/(27r), F reads 

^a(JV. JV) ^ ^ (2^^ + 2]v)) 

N,N fL,rh ^ ^ 

(2.3) 

Here a{N,N) = h^hfj where the are well known natural numbers. 
They represent the exponentially increasing degeneracy of the oscillations of 
the string and are the origin of the Hagedorn temperature in the canonical 
single string partition function q{j3). d{x) is the Heaviside function whose 
value is 1 if X > 0, zero if x < 0 and 1/2 when x = 0. As expressed in (j2.d|l . it 
is easy to see that when R is big, T{E, R, 1) converts into a multiple sum in 
which only the m = 0 term (null winding) has a non vanishing contribution 
because the Heaviside function in (Q vanishes (for m 7 ^ 0) when E < R/a'. 
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So, we have 


r(B,fl.l) = ^ e (e'‘ - --m) {E<Rla') (2,4) 

Nn^ a / 

This can be seen as T compnted for an inhnite nnmber of relativistic 
particles with squared masses AN/a'. There is a total polarization degeneracy 
of value bff and the momenta are fully discrete. It is also noticeable that 
the left-right level matching condition (a fully closed-stringy ingredient) has 
become as trivial as in the open space case. Anyway there still remains a 
stringy ingredient in the collection of particles as a whole that is reflected 
in the exponential growth of the coefficients ^at. This is the string seen as a 
collection of quantum fields or analog model. 

Now it is direct to use the standard argument to compute the sum as 
an integral based on the fact that counting the number of accessible states 
corresponds geometrically to compute the number of points of a cubic lattice 
inside a sphere defined by the argument of the Heaviside function in (ED). 
This is just accomplished by integration. This counting approximates well 
the sum as long as the radius of the sphere is big enough to contain a big 
number of points so as to minimize the relative error that appears because 
of the clear bad counting near the spherical surface. We will see in section 
4 a precise way of computing this error. The result of counting points by 
making an integral will finally give a term proportional to the volume V. 

It seems common sense to apply the same reasoning to the crude r{E,R,l). 
However, there appear difficulties. Phase space geometry is slightly more 
complex but, at first sight, tractable as the volume enclosed by a kind of 
an elliptical deformation of an eighteen-dimensional sphere (to be concrete, 
we focus on the superstring) producing two groups of nine semi-axes each of 
squared length + 2N))R‘^ and the T-duality transformed (i.e, 

what is obtained under R —> a'/ R) respectively. What spoils the reasoning 
of counting states by computing volumes by integration is the closed string 
level matching condition that makes the lattice one that cannot be built upon 
replicating a cubic cell. However, with respect to the problem of knowing 
how the volume dependence works in relation with computing integrals, the 
fictitious system obtained by putting by hand N = N instead of the con¬ 
straint in (Q will serve. The reasoning on counting states by integrating 
to get a volume is then immediate and we may try to forget for a moment 
what the error linked to surface effects could be, although we anticipate it 
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will be crucial. It is clearly easy to realize that the phase-space volume of the 
elliptical object with the given axes is independent of R because, just with 
a simple change of variables with unit Jacobian, we make R disappear from 
the argument of the Heaviside function. So, surprisingly enough, when the 
integral is a good approximation to the sum, no volume dependence appears. 
Now, if we put this result and the computation of the big volume limit (which 
gives the same result for the closed string with constraint and the hctitious 
trivialized system) together, it is clear that the difference between the sum 
and the integral, i.e. what could be called the rest from a point of view 
approximating sums by integrals, must be where volume dependence lies. It 
is then mandatory to study the exact relation between computing sums by 
means of integrals and sums themselves. This will be equivalent for us to 
studying the Euler-Maclaurin summation formula in a multivariate form. 

3 The Euler-Maclaurin formula to make sums 
by means of integrals 

Euler-Maclaurin summation formula was conceived and has mainly been used 
as a celebrated formula for numerical integration by converting integrals into 
sums. This formula was discovered independently by Euler and Maclaurin 
around 1740. It can immediately be seen to be an extension of the trapezoidal 
rule for integral approximation of one single variable functions. We need it 
here to make the converse interpretation; i.e., to compute sums by calculating 
integrals. Consequently, we need a multivariate generalization. 

Even having found a reference to the existence of such a multivariate 
generalization in |3], we have not been able to hnd the equation giving it. 
So, in the appendix, we present a generalization based on a systematic ap¬ 
plication of the Gauss theorem to a natural extension of the beautiful proof 
for this formula presented in j^l (see also [HI and the appendix here) based on 
the periodic Bernoulli functions of one single variable. The main idea roots 
on looking at the problem in a distributional sense^ and using the periodic 
Dirac delta function as the natural mediator to convert sums into integrals 
and vice versa. The Euler-Maclaurin summation formula will emerge after 
realizing that the periodic Dirac delta function in n-dimensions will produce 

^Although it can also be understood only with the need of the Stieltjes integration as 

inig. 
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one of the distributional generalizations of the Bernoulli periodic functions. 
Being more explicit, 


5per(^) = ^ 5{x-k) (3.1) 

k&Z" 

in such way that 

J dxg (x) {f (x)) 6per{x) = ^ j dxg {x) d {f {x)) 5{x - k) 

= ^'gik)= g{k)d(j\k)^ (3.2) 

k£D k£Z^ 

where ^ is a summation modified by taking only half oi g{k) when k G dD, 
the surface dehned by the equation f[x) =0 which bounds the compact 
region H as a subset of M”. The one half factor comes from the dehnition of 
the Heaviside function of one variable given above. It is worth to stress that, 
when a product of more than one Heaviside function is necessary to dehne 
the domain D, the prime will mean a certain combination of one half terms 
for lattice points on the boundary (see the appendix). Now things start to 
produce a multivariate Euler-Maclaurin formula because 

1 - Bq{x) = (5per(T) (3.3) 

Here Bq{x) is obtained from Bi{x) from the relation 

Bo{x) = V-B,{x) (3.4) 

where 


k = -— V ^ 

27r ^ 

fcG(Z")* 

which is an expression of Bi{x) thorough its Fourier expansion that empha¬ 
sizes its real character. 


(3,5) 

k^ 


Bi(x) = 


2m 


E 


^2'Kik-x 


k^ 
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An exact formula can now be written in which it is explicit how a multiple 
sum can be computed thorough integrals, to wit 

9 {k) = f dx g{x) - f dS ■ (^g {x) + f dx(yg)-Bi (3.6) 

keD D dD D 

Where we have used that 

I dxB, ■ (V/) 5 (/ (x)) gix) = -J dS-(^g (f) (3.7) 

R" dD 

with dD the boundary of the region enclosed by the surface dehned by / (x) = 
0 with a unit vector normal to the surface given by Cg = — (V/)/|V/|. By 
simply taking Bi = —V//|V f\ and g{x) = 1, (I3.7|l is also useful to get Ids'! 
when the volume element is known because it is valid for any Bi and g{x). 

In the appendix we will show how to obtain more terms by using the other 
multivariate generalizations of the Bernoulli functions given by i?n(even) and 

Bn (odd) • 

4 Simple applications to get r(^, i?, 1) for par¬ 
ticles and strings 

It seems justihed to start this section satisfying our curiosity by the simple 
application of Euler-Maclaurin’s formula to the relativistic classical gas. We 
suppose we have a single particle (of mass M) with accessible energy levels 
E‘^ = M'^ + /R"^ in d spatial dimensions. The lattice of momenta is the 

one generated by the integers, i.e. The equation dehning the domain D 

can be taken as /(x) = {E^ — — x^ = 0 where x = |x|. It gives a 

sphere in phase space of radius p = R{E‘^ — The vector Cg is the 


^If, after several applications of the Gauss theorem, the integral over D representing 
the rest is suppressed, one finds that the remaining terms generate an asymptotic series. 



unit radial vector. With all this we have 


27r‘^/2 


^ p-'+ (27rkp) 

n k^O 


p ->+oo 


27r‘^/2 

d T{d/2) 


+ ^ fc-(<=*+i)/2 cosUnkp- 

MO ^ ^ ' ' 


(4.1) 


We see that, as expected, the absolute error coming from the surface term 
grows with the radius of the sphere, i.e. grows with energy and R, but the rel¬ 
ative error of the surface contribution to the phase space volume of the sphere 
decreases with energy and R. In terms of the conhguration space volume V, 
we see that, when p is big, the relative error goes as \j\fv. 

We remark that the phase space surface term (what would be the error when 
approximating sums by integrals) oscillates changing its sign. It is important 
to stress that, for a single object with several discrete mass levels, whenever 
a mass channel opens, there is a relevant contribution from the surface cor¬ 
rections because p < 1 around that mass. The contribution just vanishes 
when the channel opens. 

The next application will be to the computation of the number of acces¬ 
sible states for a single string for the hctitious trivialized system of section 2 
with N = N and unrestricted winding and momentum vectors. In this case 
g{x) = 1 and consequently Vp = 0. The domain D is the region bounded by 
the surface given by 


^ -AN/a')^ R^{E^-AN/a') 

Here we keep N hxed because we are interested in studying the dependence 
on the volume V by summing hrst on windings and momenta. Since g{x) is 
constant, the multivariate Euler-Maclaurin formula just gives the difference 
between the sum and the integral in terms of the hrst Bernoulli function as 


^ — dm dn = — dS ■ Bi 

(rn,n)£D jy gjy 


(4.3) 
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The integral over D is the volume of the ellipsoidal region in the 18- 
dimensional phase space of momenta and windings that, as stressed in the 
introduction, is clearly independent of the radius of compactihcation. 

Our challenge will be to show (rather to prove out) that, when taking 
the limit of big volume, the surface term, that gives the error made when 
substituting the sum by the integral, is able to produce a contribution that 
cancels the integral approximation and give also the decompactihcation limit 
term that goes as V. We will take d = 9 and focus on the superstring case. 
After all, the bosonic string can only be used in this context as a sort of a 
toy model. 

To write the surface integral, we take polar spherical coordinates in both 
the nine dimensional space of the variable rh and the one for the variable n. 
Next, we take n = \n\ = Rr cos 6] m =\m\= r{a'/R) sin 6 ' with 6 G [ 0 , 7 r/ 2 ]. 
This turns dOl) for the ellipsoidal surface into the equation for a sphere of 
radius — AN/a'Y^"^. The surface element is given by 


d^ 


d6* dhlm dhlr 


= a 


/9 


- 


4iV\ 


17/2 


a' J 


1 sin® 0 cos® 6 


R"^ . 2 n cos^ 6* ^ 

— sm 6 + —^ Cs 
a ^ R^ 


where ffn,m are the solid angles in the momentum and winding variables 
respectively. The normal to the surface unit vector is given by 


- 


R^ .2/1 cos ^ 9 ' 

sm^6' 


a 


12 


R‘ 


\ ^ R . „ ^ cose\ 


with Cm, n the unit radial vectors for the space of the continuous momentum 
and winding variables respectively. 

We use (j3.5|) to get dS k ■ Cg sin (27r/c ■ x)/k‘^ for a hxed mode k = 
(/Cm, ^n) and, after that, sum over /c 7 ^ 0. The surface integral equals 


27r 


(sin 6 'ni ^7 sin 


(fcm , 


X (sin^ cos^)' 


sin 9 COS^m 7 -|- cos 0 cos 7 


X sm 


U2 U2 

( ( Ot' 

in ( 27rp ( fcm^ sin 0 cos ^ 7 -|- k^ R cos 9 cos 9-a 


(4.4) 
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where p = {E‘^ —AN/a'Y/‘^, and are the solid angles for the eight¬ 

dimensional subspaces of the winding and momentum variables orthogonal 
to the g and Xn 9 axes respectively. We have chosen the hxed vector k 
along a hxed direction taking along the X^ g axis and k^ along X^ g. 

If we expand the sine function in the last factor of (D as sin(27rpfcm^...) 
cos (27rpfcn i?...) -f cos(..) sin(...) we see that when fcn 7 ^ 0 and pa' jR -C 1 
the sinus function of the 1/R term varies slowly as long as the cosine does 
rapidly producing a small contribution to the integral. So, dually, in this 
limit (formally, R —> -|-cxo) only the fcn = 0 contributes signihcantly and we 
can write for the integral 


—a'^p^^ [ d^ 7 (sin6'm 7sin6'n 7 )^ sin®6' 

2 77 

ktn^O 


8„COs6'm,7 . 

X COS o — - -sm 

km V 


(271 p kj^ 


R 


sin 9 cos 9 


m,7 


(4.5) 


The next step is to realize that whenever pa'/ i? 1 the sum over km 7 ^ 0 
can be substituted by an integral. Since the sum does not include the nine 
dimensional zero vector, we should exclude the region in the lattice around 
the origin bounding it by a cube centered in the origin. It is worth to note that 
one must be careful with the points of the lattice that belong to the boundary 
of the cube because their contribution will be corrected by powers of one half 
factors coming by the vanishing arguments of the Heaviside functions dehning 
the domain inside the whole space. Schematically what we have is 


V- / 

(R9—cube) 


-|- corrections for points that belong to the boundary 


(4.6) 

In other words, the integral would approximate a sum of the primed type 
previously described. To compute a sum without any restriction over the 
points just on the boundary, we must correct the integral with the adequate 
contribution for these points (see the appendix). The integral taking out 
the cube can be written as the difference between the integral to the whole 
space minus the integral over the cubic region. This trick allows us to take 
advantage of the spherical coordinates to compute for the whole space and 
the Cartesian coordinates when calculating over the cube. 
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By making the sum over k-^ as an integral including 0 just after the other 
integrals are performed, we get that (0^ equals the volume of a nine dimen¬ 
sional sphere of radius i.e. 27r®/^(i?p)®/(9r(9/2)). Now one has to com¬ 
pute an integral over the region bounded by the cubic surface centered in the 
origin. In this compact region, |fcm| < v^- This together with the fact that 
R is big makes sin (27rp fcm a' sin 6 cos 6 -^j/R) ~ 271 pk^ a' sin 9 cos Omj/R a. 
good approximation as long as R formally goes to inhnity. 

The integral over the cubic region can then easily be done using Cartesian 
coordinates. We have 


+ 1 +1 

a'V® [ d/cm,!--- [ dfcm,9 [ d9 d9,n ,7 d9n ,7 {sm 9^ j sm 9n , 7 ) 


-1 


X sin^° 0 cos® 6 * cos"' 6 *ni 7 = 


i8r(9) 


(4.7) 


The presence of a' means that this term is of stringy nature. The factor 2® 
is the only contribution coming for the integrals partially representing the 
sum over k^- Now we need to correct this result by adding the contribution 
of the points on the boundary. The rule is clear (see the appendix): for the 
lattice vectors with a single component equal to ±1 we have to add a term 
l/2(/(l, 0,..., 0)-|-/(—1, 0,..., 0)-|-permutations), for the lattice vectors on the 
boundary which have two components equal to ±1 we have to add other term 
l/2^(/(l, 1, 0,..., 0)-f/(—1,1, 0,..., 0)-l-/(—1, —1, 0,..., 0)-|-permutations) and 


so on. The function / is in fact a constant because inside and over the 
cube the sinus function can be approximated by its argument when R is 
big enough. For each vector we then have a constant contribution given by 
27r® q;'®P^®/( 18r(9)). The combinatorics involved is not too difficult. The 
vectors that contribute 1/2 of the constant value are the ones with null 

f9\ 

components but one single -|-1 or —1; there are 2 W 1 such vectors. The 
ones which give 1/2^ of the constant value are the ones with two components 


equal to ±1; there are 2^ 



of them, and so on. The total correction from 


12 



all the lattice points on the cubic boundary is 


T8m“ gy 


27rV^® 

18r(9) 


a 


/ 9 


( 2 » 


1 ) 


(4.8) 


Consequently, we have shown the way the volume reappears if one starts 
from a regime in which the integral is a good approximation to the sum 
to compute the number of accessible states for a single object. The whole 
calculus for the trivialized system can be summarized as 


>^■(9 ( - -m^ - -- 


n , ^2 R^ ^2 1-2 4iV 

dmdn -diE - 77:m —-- 


a 


/2 




a 


Ellip. 


pa'/R<^1 


Q,/9pi8^9 ^ 32 (pR)^ 


362880 


+ 


945 


362880 


dS-Bi 

( 2 ‘ 


29 /o9 


1 )) 

(4.9) 


We are now prepared to confront the pure string system. The complica¬ 
tion added to the trivialized system is the left-right level matching condition 
that is a function of winding and momentum numbers. We have to hnd a 
representation of the Kronecker delta in ( 1221 ) in a way that can be used by 
the Euler-Maclaurin formula thorough a function of real variable as g{x) in 
(IS 2 D. To this purpose, we represent the Kronecker delta as an integral 


+1/2 

/ , / ^ . /A - -NX r sinlvT (A — m • n)) 

dt exp (-2 7rit (A - m ■ n)) = SA,rn-n = - tt -^— (4-10) 

TT [/a — m ■ n) 

- 1/2 

which holds for A = N — N integer and rh, n integer vectors. 

The number of accessible states for one string which has energy lower 
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than a fixed value E can then be written as 


+ 1/2 

=^a(iV,iV) ^ j dt 

N,N n,rn 

For given oscillator numbers, we can represent the sums over windings 
and momenta as integrals by using the multivariate Euler-Maclaurin formula 
to get 


+1/2 

- 1/2 Fm 


a' 

+ 1/2 


- ^ (4iV + 2A) ) 
a' J 


= a'^ I dt e 


—2 Tvit A 


- 1/2 


dmdn e 


2 TV it a' E •n 


-2 


dm dn e 


D 


2 7Tita' E‘^ m-n 


Rm n\ -* fEa'm 


■ Bi —— , ERn 


a' ' Rj \ R 


X E ^5(1 — m^ — — 


2 ^2 


4iV + 2A 


E‘^a' 


+ 27iitE / dm dn e^ 


Rn, ) ■ Bi 


Ea'm 

R 


, E Rn 


D 


(4.12) 


where the domain D is the region bounded by the surface defined by the 
vanishing of the argument of the Dirac delta that appears because it seems 
convenient to represent the surface integral as a volume integral over the 
whole (continuous) phase space. 

Now one has to compute the limit of F when the volume gets big to see 
the way it finally becomes proportional to the volume of the hypertorus as 


14 



one departs from a situation in which no volume dependence appears because 
the first integral is a good approximation. When energy is high enough, one 
gets that D, the domain of integration^, is such that the number of points 
inside the boundary is really big when compared to the error over the surface 
that increases, but with a lower power of p'. 

Let us now give the detailed proof of the way the high volume limit shows 
up from (ICTl) . The hrst step is just to expand as a series of 

powers of its argument. The zero order in this expansion gives, for the hrst 
and second integrals, the same contribution as for the hctitious trivialized 
toy model (after changing p by p'), so the analysis made for this example 
applies for it and hnds here its real purpose. The contributions from the 
odd powers vanish^. For the last integral in (iron . it is easy to see that the 
contributions from the even powers of the expansion of the exponential also 
vanish. 

We have to compute a sum over in the integrals containing Bi. This 
sum can be computed as an integral over the real variable k plus the con¬ 
tribution from taking out the point fcm = 0 (see ( 021 )). The contribution 
integrated over /c G coming from the (2j — l)-th power in the third in¬ 
tegral exactly cancels the part integrated over k from the 2j-th power of the 
second integral (always in the limit in which p a'/R 1). 

Finally, for both sums over fcm, we are left with the contributions from the 
points on the boundary of the cube centered at the origin plus the integral 
over this domain (cf. (jd.bjl L It is easy to see that these contributions from 
the 2j — 1-th order in the third integral and those from the 2j-th order in 
the second one add up to exactly cancel the 2j-th order of the first integral, 
which is the part manifestly independent of the configuration volume V. This 
is the rigorous proof of the connection between the particle gas regime and 
the volume independent situation for F(F^, R, 1). 

'^D is always an ellipsoid with semi-axes given, after putting p' = {E'^ — {AN -|- 
2A)/a;')^^^j by p'a'/R and p'R 

®We could have written cos 2 tt t (A — m • n) instead of the exponential as the integrand 
in to make this manifest. 
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5 The limit of decompactification for the gas 
of strings 

Now it is mandatory to try to get the first consequence of what we have 
learned. We think that the issue on whether the uncompactihed regime for 
the gas of strings can be recovered from the compact description deserves 
a certain reanalysis after putting together what we know now on T{E, R, 1) 
and what is presented in [3]. 

By using the inverse Laplace transform techniques the best that has been 
obtained for the density of states of the gas of strings with compact dimen¬ 
sions is a series that converges whenever [2] 

E > = Eh (5.1) 

{d is the number of big dimensions). The term on the right hand side of this 
inequality is a sort of Hagedorn characteristic energy. The important point 
here is that in j2] it is explained that, when R is big enough as compared with 
a/^, K can actually be computed to give a precise physical dehnition of the 
Hagedorn energy Eh as the energy at which the system reaches the Hagedorn 
temperature or approaches it from a low energy regime dominated by a gas of 
massless objects. Hence, the Hagedorn regime will be dehned as the regime 
at which the massive modes appear in the gas. These modes have masses 
that are proportional to l/\f^ which is the characteristic mass scale of the 
stringy effects. They are then oscillators, whose mass is a positive integer 
times 1 / a/^, or windings of mass R/a'. To be more concrete, when one has a 
system with an Avogadro’s number of strings, looking at its Hagedorn regime 
dominated by oscillators means looking at the system when the the energy 
is such that E/Nh > 2 j\/c/, where Nh is the number of strings that the 
gas presents at the Hagedorn energy as given by the characteristic Hagedorn 
temperature (cf. jS]). 


Eh = C,oVP],^^ (5.2) 

where Ciq is a known constant. 

The bound 2/\f^ is precisely the minimum mass of a non vanishing 
vibrational mode of a single string. At this point, we must take into account 
that there is a relationship between the number of strings and energy and 
also volume (or R). It is the result of the fact that our microcanonical 
description is one very special (it is not rigorously microcanonical) for which 
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the chemical potential /i vanishes. This condition for equilibrium gives in 
that Nh ~ which simply states that the number of strings at this 

characteristic energy is hnally proportional to the volume, like in the black- 
body in equilibrium with a gas of photons. This is so because one can always 
think of the system as reaching or approaching the Hagedorn temperature 
Th at the energy Eh departing from the low energy regime in which strings 
are in their massless modes. Thus, Eh is the energy which separates the low 
energy regime of the gas dominated by the massless modes of the string (the 
a' —0 limit taken on the mass formula) from the string or high energy 
regime in which vibrational or winding modes appear. If R is big enough, 
there is always a low energy regime with massless modes for the gas so we 
now understand why the condition E — Eh > 0 appears the same for different 
ab initio topologies. This was seen as something notorious in the past. 

If one naively thinks that the condition to obtain the uncompactihed 
result is to impose E < R/a', i.e., the replica of the condition over the 
single string, but now E stands for the total energy, it is easy to see that 
being in the Hagedorn regime together with this condition forces R to be 
small instead of big as it was supposed. This way, it would seem that there 
is no domain of energy and radius for which the characteristic asymptotic 
behavior in uncompactihed space can be recovered. The physical picture of 
this behavior at the Hagedorn regime is that of a single string absorbing 
all the energy in its vibrational modes and approaching equilibrium with 
a sea of low energy strings at the microcanonical Hagedorn temperature. 
It is precisely the physical picture what makes clear that, to obtain the 
Hagedorn regime dominated by the vibrational modes characteristic of the 
uncompactihed result, the condition to add to that of being in the Hagedorn 
regime [E > Eh) is E — Eh < R/a'. It simply states that the single fat string 
cannot have winding modes. Now, this is compatible with the fact that R ^ 
V^. In other words, to get the fat string dominated regime characteristic of 
having open spatial dimensions, one needs a radius big enough so as to have 
that the energy of the system minus the Hagedorn energy does not suffice 
to get windings in the gas in equilibrium (or quasi-equilibrium). Vibrational 
modes are much lighter and, after reaching Eh, massive strings appear to 
stay in meta-equilibrium with the sea of massless modes. It is worth to 
notice that the volume can be very big but the momentum that a string in 
the gas can have still belongs to a discrete set; we then have and effective 
decompactihcation but not a pure one that would be the one in which the 
nine torus certainly becomes and the momenta get dense. We will treat 
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this point later when commenting on the Jeans instability problem. 

What happens with the approximate series computed in j2] is that it is 
adapted to describe the Hagedorn regime when E — Eh > R/a'. Mathemat¬ 
ically this is analogous to the way in which the integral approximates the 
sum in F to only describe the high energy approximation. 

If, from the very beginning, one has that R is of the order of a/o^, then the 
mass formula tells us that there is no radiation dominated phase (the stan¬ 
dard black-body regime with massless objects) because a massless object 
with non vanishing momentum is indistinguishable from a massless object 
with a non vanishing winding. However equation (EH) is still valid with 
d = 0 giving that high energy means E > K/\f^ (cf. the way this condition 
appears in j2]). In any case, the total energy has to be bigger than the mini¬ 
mum momentum or winding mass which is as big as the mass of a vibrational 
mode. Of course, if winding modes are present in the gas, from the expression 
approximating Q{E,R), it seems that one cannot get the characteristic neg¬ 
ative specihc heat phase for the gas with uncompactihed dimensions. But, 
contrary to what is stated in j2] this does not in any way imply that there 
is no domain in the plane {E, R) for which the uncompactihed behavior can 
be found. 


6 Comments and outlook 

Before commenting on the contains of the work, may be some more general 
remarks come into place to justify the interest of the topic we have dealt 
with. First, it is important to notice that perhaps one of the most active 
research topics in String Theory has been the one about the Hagedorn phase 
which is characteristic of a gas of free and perturhatively interaeting strings. 
Several published works have claimed to have completely solved the problem 
mainly in two ways: there is really a (may be hrst order) phase transition 
at the Hagedorn temperature or the Hagedorn temperature is a maximum 
one. Things are more subtle than some statements that have been published 
and are being taken as true. This is the motivation for us to pursue a plan 
to unambiguously settle what is really well established and what is simply 
wrong when treating the topic. To this purpose we have already clarified 
the way quantum statistics effects are taken into account in the problem, the 
connection between the presence of non extensive terms and negative micro- 
canonical specific heats, and the physical meaning and value of such peculiar 
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systems jS]. Systems having a negative microcanonical specific heat are a 
subject of intense active research by the community of theoreticians working 
on Statistical Mechanics and by atomic and condensed matter experimental¬ 
ists. 

The main criticism that can be made to any work on the free string gas is 
very well known because one can find it, for example, in a classic review on 
String Theory as it is the two volume work by M. Green, J. Schwarz and E. 
Witten [7] . It is a criticism that simply brushes away the Hagedorn problem 
because Jeans instability and the nucleation of black-holes simply make ill- 
defined the canonical treatment when interactions are plugged in. It is then 
often argued, when treating thermodynamical properties of string gases, that 
no thermodynamic limit can be taken because we are dealing with systems 
that include gravity, i.e. long range interactions, and no infinite volume 
(thermodynamic) limit exists (although there is certain controversy on this 
subject, it seems that the theorem presented in m is an exact proof of this 
fact for classical gravitational interactions). As we are unable to take this 
limit, the concept of temperature and the idea of a possible phase transition 
or having a maximum temperature must be discarded because we will be 
unable to follow the ideas coming from Statistical Physics (see again [7j). 

This statement is both wrong and true; it needs clarification. When no 
thermodynamic limit exists, the canonical description, whose main param¬ 
eter is temperature (T), becomes ill-defined. The reason is that, from the 
mechanical point of view, T is a derived parameter, which has no direct 
significance. On the other hand, a microcanonical description relies upon 
the use of energy E, a quantity conserved for an isolated system and well 
defined in phase space. Entropy has also a clear meaning as the logarithm 
of a volume in phase space (we do not really need a probabilistic definition: 
S = klnW versus S = Pd^Pi)- Statistical Thermodynamics, extensiv- 
ity is assumed and the thermodynamical limit (as V oo) is taken. It is 
in this situation that both canonical and microcanonical descriptions appear 
to be equivalent. A phase transition will show as a non homogeneity of the 
systems, a non extensivity. 

In our problem, even when the string coupling vanishes {qs = 0) and we 
do not have long range interactions, we do not have an extensive entropy 
[ 5 . It can also be shown that microcanonical and canonical descriptions are, 
in some cases, inequivalent. After introducing interactions in our system, 
no thermodynamic limit can be taken as the result of the anti-screening 
properties of gravitational forces; we will then be forced to treat the system 
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microcanonically, this being the ensemble adapted for ’’small” systems. This 
would let us also study phase transitions in hnite volume systems which, 
by the way, can also be seen experimentally in ’’small” systems. A very 
clear exposition of these facts and its relation to negative specihc heats can 
be found in jH] where, among others, classical self-gravitating systems are 
specihcally treated. 

For the particular case of self-gravitating systems, where Jeans instability 
appears, an important effort has been made to treat the problem. In the work 
by Gross-Perry-Yaffe jH], where an approach from the quantum held theory 
point of view was presented, the need for a microcanonical treatment was 
remarked. Since then, the traditional reasoning to avoid Jeans’ instability is 
to put the system in a hnite volume as we have done in our work by using 
a hypertorus. Then the procedure is to rapidly scale the string coupling 
constant to zero keeping the canonical energy density constant so as to allow 
the volume to become large. Alternatively, one can hx the coupling as very 
small and restrict the volume to be large but hnite®. It is standard lore that 
to solve infrared divergences due to long range interactions one can use a 
hnite volume. With hnite volume there is no = 0 momentum for the 
massless objects. Jeans’ instability really appears thorough the computation 
of a non vanishing contribution for the temporal components of the vacuum 
polarization tensor with p® = 0 and p —0 in the pure decompactihcation 
limit in which the momenta become dense, in particular, around zero. To 
avoid Jeans’ instability in the canonical description, one can have a really big 
radius (volume) as compared to a/Y, but the momenta must still be discrete. 
This is the precise meaning of ’big but hnite’ and the diherence between an 
ehective decompactihcation (big R) and a pure decompactihcation in which 
the torus becomes and the momenta become then dense. 

A very remarkable fact is that, with an up to now complete ignorance from 
the string theorist community, the treatment of self-gravitating systems has 
recently been a very active research held in Statistical Physics jT^j. Although, 
in these works, gravity is treated classically (as Jeans did), very interesting 
results appear related with phase transitions for the system stemming from 
negative specihc heats and the need for a short distance cut-oh (which, by 
the way, naturally appears in String Theory and should be provided by it). 
These results are largely seen to be independent of the way the cut-oh is 

®More concretely one has that E < ^ R‘^~‘^/Gn as the condition 

to avoid Jeans’ instability 
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introduced na. For a review see in 

We believe these works signal the importance of a microcanonical treat¬ 
ment for string gases even in the free case. In this work, we have set aside 
questions related to Jeans instability because we hnd it compulsory to make 
a microcanonical treatment of the problem from hrst principles in String 
Theory paying special attention to its relationship with the work in ca 
(and references there in). In a completely independent way, interesting work 
na has been made on the constraints that the classical Jeans instability (as 
it is gotten from a quantum held theory in the canonical description) and 
black-hole nucleation have upon the possibility of really having a Hagedorn 
’’visible” phase that can actually produce black-holes from highly excited 
(fat) strings at higher energy. However, no direct calculation in the micro- 
canonical description of the Jeans effect has been done. Knowing as we do 
the results recently found about the microcanonical description of classical 
self-gravitating systems, we, by no means, hnd ourselves justihed to simply 
apply the results about Jeans’ instability gotten from a hxed temperature 
description. May be we are being too much conservative, but there are too 
many subtleties in the subject that deserve some precaution. Finally, to un¬ 
derstand how Jeans instability appears related to the string loop expansion 
computation of the canonical free energy, the work in ng includes a expla¬ 
nation of the relationship between the perturbative loop expansion and the 
appearance on Jeans’ instability in the computation in [^. 

There can be another criticism about the intrinsic instability of the toric 
box seen in a cosmological context. Problems like decompactifying dimen¬ 
sions, the running dilaton, etc. that are characteristic of toric zero temper¬ 
ature vacuum solutions have to be posed, at hnite temperature, in a cosmo¬ 
logical picture where the stress tensor of the string gas enters in the Universe 
evolution equations like it was done long time ago in uni- Whether a torus is 
a good background can then only be seen afterwards although lots of works 
have been published and are actively being produced in this direction (see, 
as a single example, ini). Anyhow our main interest in this work is pre¬ 
cisely decompactihcation, although non-dynamical. Sections 4 and 5 have 
been devoted to show how the toric box can decompactify and we can reach 
a situation where, for the string gas, all the dimensions are effectively open. 

Let us now focus on what has been the precise subject of this article. By 
computing F(U,i?, 1) for a string in compact space (or its derivative with 
respect to the energy) by applying a high energy limit analogous to the one 
in the regular particle gas, one simply hnds a limit in which the sum upon 
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windings and momenta can be converted to an integral with big accuracy. 
This integral is easily seen to be independent of the volume. After all, con¬ 
verting sums into integrals is a way of getting a more explicit dependence 
on the volume than the one seen just in the raw sum. Asking why not to 
use the crude sum instead of the integral has to do with knowing the volume 
dependence for the thermodynamical system as it happens in the classical 
free particle gas. What is dropped in this approximation, the phase space 
boundary contribution, just contains the dependence over the volume that is 
not present in the high energy limit. The approximation is only valid when 
energy gets big and the radius of the hypertorus is relatively kept not too big. 
When the radius grows and so the relation a' E/R is small, the boundary 
contribution increases its weight and cannot be dropped, it actually produces 
a term proportional to the volume as in the open space case as long as R 
becomes bigger with hxed energy. In other words, studying hnite size effects 
in a string system is a subtle issue because of T-duality. 

To better see this, we have shown that what is easily seen true as ex¬ 
pressed in a set of variables is also true when things are re-expressed in a 
different set of variables, whenever all the contributions are taken into ac¬ 
count. In the problem at hand, this statement means that the high volume 
limit so effortless computed over r{E,R,l) as written in (I2.d|l can be re¬ 
covered when the sum over discrete momenta and windings is represented 
thorough integrals as though they were continuous variables. Furthermore, 
as long as we keep using the representation of the periodic delta function 
in terms of phase space coordinates without introducing more Bernoulli pe¬ 
riodic functions, the task of getting the high volume limit is still a simple 
one because the integral representing T will hnally get a contribution from 
only one of the Dirac delta functions in the dehnition of 6 pej.{n,rh). Namely, 

in the limit R —>• -|-cxo of the product'd (^p' ^ ^ j 5per(ff, hi) only 

the contribution for 6 {rh) survives. It is now obvious how the high energy 
calculations giving an independent of the volume r{E,R,l) are related to 
substituting the sum over windings and momenta in the computation by a 
multiple integral. What is now new is the fact that the volume independent 
regime can be connected in a smooth way with the regime proportional to 
the volume when the equality between the sum and the integrals expressed 
thorough the EM formula is taken into account. Furthermore, there might 
be intermediate regimes as for the volume dependence is concerned. As an 
example, if there really is a regime proportional to the area of the hypertorus 
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and it is compatible with a Hagedorn regime in which one single string dom¬ 
inates the system with nnll chemical potential, then the many-string system 
wonld have an entropy with a term proportional to the area of the snrface 
enclosing the volnme V. 

Going over to the technical side, there are certain ambignities when rep¬ 
resenting a snm by integrals. The hrst one is related to the election of the do¬ 
main D which contains the lattice points we are snmming over. For example, 
it holds that J^kez 'd(/+0.8 —fc)'d(fc — (n— 0.2))'d(/—n) = 1 = 

{Bi{l + 0.8) — Bi{n — 0.2)) which certainly eqnals / — n -|- 1 (to simplify, n , I 
are natnral nnmbers, {1.8) = 0.8) and it also eqnals J2k&[n-ox i+oy] This 
shows that there are, in general, more than one integration domain and that 
issne affects the valne of the integral over D as the hrst approximation to the 
snm. However, this ambignity is harmless becanse it rehects the tolerance 
from the distance of a point in the lattice of integer components to the nearest 
neighbor. This aspect is negligible in the calculation of T. What seems more 
relevant at hrst sight is the ambiguity in representing the Kronecker delta 
function expressing the closed string constraint. It is immediate to see that 

sin 7r(A — m ■ n 

7 T {A — rh ■ n) 

is another representation of 5 a, mn for m, n integer vectors that could be 
used instead of the one in (ICTl) . It is manifest that, because of the cosine 
term and after applying the Euler-Maclaurin summation formula, they are 
diherent as real variable functions although coincide for integer arguments. 
Taking the limit a' E'^ —> -f-cxo over the hrst integral upon D on the right 
hand side of for both elections is a good way of estimating the ehect 

of changing the representation. We get that 


the integral dt cos 27rt(A — m ■ n) 


cos 7 r(A — rh ■ n) 


+ 1/2 

lim a' E'^ / dt cos2Trt(A — a' E'^ m ■ n) = 5 i - - — rh ■ n] 

J ^ > \a'E^ J 


- 1/2 

1 


lim a'E'^ [ dt cos2TTt(A — a'E"^ m ■ h) = - 5 ( - - — rh ■ n] 

J ^ ’2 \a'E^ J 


where we understand that, when A 7 ^ 0, we will have contributions from 
terms with A of the order of a' E^. It is clear that after summing up the 
contributions from the three integrals in (sm, the result of computing the 
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sum by making integrals is independent of the way we represent the Kro- 
necker delta. 


Note added 

After disclosing our work, we have found ^H] (and references therein) in which 
several multivariate generalizations of the Euler-Maclaurin formula are pre¬ 
sented stressing their usefulness to numerically compute integrals by sums. 
It seems clear that the restrictions this numerical problem imposes makes 
this project a more difficult mathematical task than it was ours. 
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Appendix 

Let us present the main ingredients of the Euler-Maclaurin formula as ob¬ 
tained in [S] (see also [B]). 

Bernoulli numbers Bn play an important role in the Euler-Maclaurin for¬ 
mula. Bernoulli numbers are concrete values of the Bernoulli polynomials 
{Bn = B*{0)) that can be dehned by the recursive formulas: 


b;{x) 

= 1 

(A.l) 

Bi(x) 

= nBl_^{x) 

(A.2) 


= 0 

(A.3) 


0 
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with n > 1 

The Bernoulli periodic functions are defined as Bn{x) = i?*((x)), with 
(x) = X — [x] = fractional part of x. Now M. Lerma in [S] realizes that the 
odd numbered Bernoulli numbers are zero except Bi = —1/2 and that for all 
Bernoulli periodic functions dx Bn{x) = 0 except for n = 0 that the inte¬ 
gral gives 1. To make things easier, it is adequate to use a modihed version of 
the periodic Bernoulli functions that has sense if Bq is a distribution instead 
of a function, but that makes the proof of the Euler-Maclaurin summation 
formula more natural and, what is important to us, easily generalizable to 
multiple integrals and sums. 

The Bernoulli periodic distributions are dehned thorough the recursive 
relations (n > 1) 


Bo{x) = 1 - Sperix) (A.4) 

B'^{x) = nBn-iix) (A.5) 

1 

J dx Bn{x) = 0 (A.6) 

0 

Now Bq{x) also verihes ()A.6|1 . The differences with the standard definitions 
of the Bernoulli functions are in Bq{x) and Bi{x). Bi{x) is now the function 


Bi{x) 



if X ^ Z 

if X e z 


(A.7) 


This can be checked in relation with (HaJ by computing Hi (x) = Jo dp -Boil/)- 
It gives 


X 

Bi(x) = x f dy S{y-k) 

fcez i 

= x — ^ [^{x — k) — 'd(—fc)] 
fcez 



if X ^ Z 
if X e Z 


(A.8) 


where it is worth to remark that the Dirac delta has been integrated using 
5 (x) = ^ 'd(a;) with 'd(x) the Heaviside function of value 1 for a; = 0. 
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The Fourier series representation of the Bernoulli functions will be ex¬ 
ceedingly useful when generalizing Lerma’s proof to multivariate sums. We 
have 


Bn{x) 



+ 00 

E 


k=—oo 

k^O 


^ 2'Kikx 


(A.9) 


The proof of the EM summation formula is now very easy and starts with 
la x) — 'Yl,a<k<b /(^) ~ /a ^xBq{x). This is not a tautology because 
we can use that B'^{x) = Bq{x) to integrate by parts on the right hand side. 
This can be done successively by using The prime over the sum results 

again from the fact that 'd(O) = 1/2. And the summation formula emerges 
beautifully. 

To generalize this to multiple sums is straightforward. We can easily 
write down the periodic Dirac delta function of n variables (see (HH)) and 
define the rest of the multivariate Bernoulli periodic functions thorough their 
Fourier expansions as two series of scalar and vectorial functions. 


n. 


(27ri 

E 


' ^2'K\k-x 

if n G (2Z) 

(A.IO) 

O’" ^ IfcF 

fcez ' ' 




P 2'K\k-x 



i^>"‘ '■ 

if n G (2Z -|- 1) 

(A.ll) 





Now 1 — Bq{x) = Yhk — and the recurrence relations are V-i? 2 n+i(a^) = 


(2n -h 1) B 2 n{x) and V i? 2 n+ 2 (ai) = {2n 2) i? 2 n+i(^) with n > 0. 

The expression we get for a sum over integer vectors belonging to a com¬ 
pact region D defined as a subset of M”, with Bi[x) and B 2 {x) explicitly 
shown, is 


9{k) = j dx g{x)- f dS-(^g (x) Bi^ + ^ j d ^-(^^2 Vfi-j f dxB 2 Ag 

keD D dD dD D 

(A.12) 
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The prime means that the lattice points /cq that are just over the surface 
boundary will contribute a factor (different from one in general) times g{kQ). 
In the case the domain D is defined as the region bounded by the surface 
defined by the subset of of points x such that f{x) = 0, a single Heaviside 
function will serve to write 'Yl,k£D ' ~ SfceZ" the prime 

will mean that the integer vectors just on the boundary surface, /cq, will 
contribute g{ko)/2 because the value of the step function is 1/2 when its 
argument gets null. 

In our applications, we have found a sum in which the domain of in¬ 
tegration is one that cannot be defined using a single step function, but a 
combination of products of several of them. More precisely, it is the sum over 
all the nine dimensional integer /c-vectors with the exception of the vector 0. 
We have chosen to take out the origin bounding it by a cube of two length 
side centered at the origin. The domain D will then be the region outside 
the cube. To get only a contribution to the sum from points in H, we must 
insert in the sum a factor 


*1^*2^...^*9 fc=l 

(A.13) 

From here it is easy to see that, for example, a vertex of the cube like 
(1,1,1,1,1,1,1,1,1) would contribute 



(A. 14) 


If one wants to compute the sum without any restriction, then one has to 
sum 1/2^ oi g ((1,1,1,1,1,1,1,1,1)) to the integrals giving the primed sum. 

In our case, we have that g{k) is constant so one can readily show that 
the integral faithfully represents the primed sum because the integral over 
the surface of the cube vanishes. 
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